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We study theoretically some symmetry properties of spin currents and spin polarizations in multi- 
terminal mesoscopic spin-orbit coupled systems. Based on a scattering wave function approach, we 
show rigorously that in the equilibrium state no finite spin polarizations can exist in a multi-terminal 
mesoscopic spin-orbit coupled system ( both in the leads and in the spin-orbit coupled region ) 
and also no finite equilibrium terminal spin currents can exist. By use of a typical two-terminal 
mesoscopic spin-orbit coupled system as the example, we show explicitly that the nonequilibrium 
terminal spin currents in a multi-terminal mesoscopic spin-orbit coupled system are non-conservative 
in general. This non-conservation of terminal spin currents is not caused by the use of an improper 
definition of spin current but is intrinsic to spin-dependent transports in mesoscopic spin-orbit 
coupled systems. We also show that the nonequilibrium lateral edge spin accumulation induced by 
a longitudinal charge current in a thin strip of finite length of a two-dimensional electronic system 
with intrinsic spin-orbit coupling may be non-antisymmetric in general, which implies that some 
cautions may need to be taken when attributing the occurrence of nonequilibrium lateral edge spin 
accumulation induced by a longitudinal charge current in such a system to an intrinsic spin Hall 
effect. 

PACS numbers: 72.25.-b, 73.23.-b, 75.47.-m 



I. INTRODUCTION 



The efficient generation of finite spin polarizations 
and/or spin-polarized currents in paramagnetic semi- 
conductors by all-electrical means is one of the princi- 
pal challenges in semiconductor spin-based electronics.^''^ 
For this purpose several interesting ideas have been pro- 
posed based on the spin-orbit ( SO ) interaction char- 
acter of electrons in some semiconductor systems ^i^i^ 
One such interesting idea is the so-called intrinsic spin 
Hall effect,^'^ ■which is the generation of a finite spin 
current perpendicular to an applied charge current in a 
paramagnetic semiconductor 'with intrinsic SO coupling. 
Such ideas have attracted much theoretical interests 
recently^"^** and substantial achievements do have been 
obtained along these lineS )^^'^^ -while at the same time 
they also raised a lot of debates and controversies i^ii^ 
A central problem related to these debates and contro- 
versies is that, -what is the correct definition of spin cur- 
rent in a system -with strong SO coupling and what is 
the actual relation between the spin current and the in- 
duced spin-accumulation in such a system,^^ In most re- 
cent studies the conventional ( standard ) definition of 
spin current ( i.e., the expectation value of the product 
of spin and velocity operators ) have been applied. As 
is well know, this conventional definition of spin current 
can describe properly spin-polarized transport in a sys- 
tem without intrinsic SO coupling. However, since spin 
is not a conserved quantity in a system with intrinsic 
SO coupling, the physical meanings of spin current cal- 
culated based on the conventional definition are some- 



what ambiguous and the actual relations between the 
spin current and the induced spin- accumulation are not 
much clear. In fact, as has been noticed in several recent 
papersr^i^i^i^ there may exist some serious problems 
with this conventional definition of spin current when 
using it to describe spin-polarized transport in a system 
with intrinsic SO coupling. In order to avoid such serious 
problems, several alternative definitions of spin current 
were proposed in these papers based on different theoret- 
ical considerations, which are significantly different from 
the conventional one and also significantly different from 
each other Another possible way to circumvent 
this problem is to study a mesoscopic SO coupled system 
attached to external leads. If no SO couplings present 
in the leads or the SO couplings in the leads are much 
weak, then the conventional definition of spin current can 
be well applied without ambiguities in the leads. Several 
recent works have adopted this strategyi^'i^ii^iii'i^ and 
some interesting results were also obtained. Of course, 
the study of spin transport in mesoscopic SO coupled 
systems is not only of theoretical interest but also might 
find some practical applications in the design of spin- 
based electronic devices<^ 

In this paper we study theoretically some interesting 
problems related to spin-dependent transports in multi- 
terminal mesoscopic SO coupled systems. We focus our 
study on the symmetry properties of equilibrium and 
nonequilibrium spin currents and spin polarizations in 
such mesoscopic structures. As is well known, symme- 
try analysis is usually of great theoretical importance 
in the study of many physical phenomena, including 
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the spin-dependent transport phenomena in SO coupled 
systems/i22, Based on the analyses of symmetry proper- 
ties of equilibrium and nonequilibrium spin currents and 
spin polarizations in multi-terminal mesoscopic SO cou- 
pled systems, some controversial issues related to spin- 
dependent transports in mesoscopic SO coupled systems 
will be investigated in some detail in this paper. Some 
symmetry properties discussed in this paper might also 
be helpful for clarifying some controversial issues encoun- 
tered in the study of spin-dependent transports in macro- 
scopic SO coupled systems. The study carried out in 
this paper is based on a scattering wave function ap- 
proach within the framework of the standard Landauer- 
Biittiker's formalism. From the theoretical points of 
view, this scattering wave function approach is in prin- 
ciple exactly equivalent to the more frequently employed 
Green's function approach in literature'^'^. The main 
merit of this scattering wave function approach is its con- 
ceptual simplicity, and due to its conceptual simplicity, 
some symmetry properties of equilibrium and nonequi- 
librium spin currents and spin polarizations in a multi- 
terminal mesoscopic SO coupled system can be more ex- 
plicitly shown. 

The paper is organized as follows: In section II we 
will first give a brief introduction of the structure con- 
sidered and the approach applied. In section III we will 
use the approach introduced in section II to investigate 
whether there can exist nonvanishing equilibrium spin 
polarizations or nonvanishing equilibrium terminal spin 
currents in a multi-terminal mesoscopic SO coupled sys- 
tem. In section IV we will study the symmetry properties 
of nonequilibrium spin polarizations and nonequilibrium 
terminal spin currents in a typical two-terminal meso- 
scopic structure with both Rashba and Dresselhaus SO 
coupling. Finally in Section V a brief summary of the 
main conclusions obtained in the paper will be given. 



II. DESCRIPTION OF THE STRUCTURE AND 
THE SCATTERING WAVE FUNCTION 
APPROACH 

We consider a general multi-terminal mesoscopic struc- 
ture as shown in FiglU where a SO coupled mesoscopic 
system is attached to several ideal leads. In a discrete 
representation, both the SO coupled region and the ideal 
leads are described by a tight-binding ( TB ) Haniilto- 
nian, and the total Hamiltonian for the entire structure 
reads: 

H = Hieads + Hsys + Hs-l- (1) 

Here Hi^ads = Spi/p and Hp 
-ipi;<p^,p^.>CT((7t^^(7p^.^ + H.C.) is the Hamiho- 

nian for an isolated lead p, with Cp-a denoting the 
annihilation operator of electrons with spin index cr at a 
lattice site Pj in lead p and tp the hopping parameter 
between two nearest-neighbored lattice sites and pj 
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FIG. 1: Schematic geometry of a multi-terminal mesoscopic 
SO coupled system. 

in the lead. We assume that the external leads are ideal 
and nonmagnetic, i.e., no any SO couplings ( or other 
kinds of spin-flip processes rather than that induced by 
the scatterings from the central SO coupled region ) 
present in the leads. In such ideal cases the standard 
definition of spin current can be well applied without 
ambiguities in the leads. Usually ( if not specified ) 
we will choose the z axis ( normal to the 2DEG plane 
) as the quantization axis of spin. Hgys = Hq + Hso 
is the Hamiltonian for the isolated SO coupled region, 
in which Hso describes the SO coupling of electrons 
and H(i = ~tT.<ri,rj>c,{Cl.^Crjc, + H.C.) describes 
the spin-independent hopping of electrons between 
nearest- neighbored lattice sites ( denoted by {vi^Vj) ) in 
the region. The discrete version of Hg^ will depend on 
the actual form of the SO interaction, e.g., for the usual 
Rashba and k-linear Dresselhaus SO coupling, one has 

(2a) 

(2b) 

where tji and tn are the Rashba and Dresselhaus SO cou- 
pling strength, respectively, ^'r^ = (C'r;,! j C'ri.j.) denotes 
the spinor annihilation operators, and A^, and Ay de- 
note the lattice vectors between two nearest-neighbored 
lattice sites along the x and y directions, respectively. 
The last term in the Hamiltonian (1) describes the cou- 
pling between the leads and the SO coupled region, 
Hg-i = -Y^ptps J2ni'^l„(7Cr„cr+H.C.), where p„ denotes 
a boundary lattice site in lead p connected directly to a 
boundary lattice site r„ in the SO coupled region and 
tps the hopping parameter between lead p and the SO 
coupled region. It should be noticed that in general the 
TB Hamiltonian will also contain an on-site energy term, 
which is not explicitly shown above. 
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Now we consider the scattering of a conduction elec- 
tron incident on a lead p by the SO coupled region. For 
conveniences, we will adopt a separate local coordinate 
frame in each lead, i.e., we will use a double coordinate 
index {xp,yp) to denote a lattice site in lead p, where 
Xp — l,2,...,cxD ( away from the border between the lead 
and the SO coupled region ) and yp = l,...,Np ( Np 
is the width of lead p ). In the local coordinate frame, 
the spatial wave function of a conduction electron inci- 
dent on lead p will be given by e~*'^^'^pXm(j/p)' where 
fcPj denotes the longitudinal wave vector and XmiUp) 
the transverse spatial wave function and m the label 
of the transverse mode. The longitudinal wave vector 
will be determined by the following dispersion relation, 
— 2tcos(fc^) + ef,^ = E, where ef^ is the eigen-energy 
of the mth transverse mode and E the energy of the 
incident electron. It should be noted that, due to the 
presence of SO coupling in the central scattering region, 
spin-flip processes ( e.g., the spin- flip reflection ) will be 
induced in the leads when a conduction electron is scat- 
tered or reflected by the central scattering region, even if 
the leads are ideal and nonmagnetic ( which is the just 
case assumed in the present paper ). Due to this fact, for 
a conduction electron incident from the mth transverse 
channel of lead p and with a given spin index a, both the 
scattering wave function \tpP"^'^(r)) in the central SO cou- 
pled region and the scattering wave function |?/'^""^(xp')) 
( Xp' = {xp' , j/p' ) ) in a lead p' will be inherently a su- 
perposition of a spin-up and a spin-down components, 

\V""^{r)) = E<'"''W4l,.'|0>, (3a) 

iv^^"'^(xpO) = EV'Sr(^p')(5I^„.'io), (3b) 

CTl 

where |0) stands for the vacuum state. The spin-resolved 
components ipaT'^ i'^p') of the scattering wave function in 
lead p' can be expressed in the following general form, 

r.rM = Spp,s,,,e-^'^'^'-^''x^M 

+ E rpZ'.'e<"'^'xi{yp'), (4) 

where (j^T^i^i stands for the scattering amplitude from 
the (ma) channel of lead p ( the incident mode ) to 
the {m'a') channel of lead p' ( the out-going mode ). 
If p' = p, the second term on the right-hand side of 
Eq.(4) will denote actually the spin-resolved reflected 
waves in lead p and (p^?^,^, denote the spin- flip { a ^ a 
) and non-spin- flip { a = a ) reflection amplitudes. The 
scattering amplitudes (f'p^i^'a' '^^'^ obtained by solving 
the Schrodinger equation for the entire structure. Since 
Eq.(|3]) is just a linear combination of all out-going modes 
with the same energy E in lead p', the Schrodinger equa- 
tion is satisfied automatically in lead p' except at those 
boundary lattice sites in lead p' connected directly to 
the SO coupled region. Due to the coupling between the 



leads and the SO coupled region, the amplitudes of the 
wave function at these boundary lattice sites ( which are 
determined by the scattering amplitudes (/)^^ ) must be 
solved simultaneously with the wave function ■(/)P'"'^(r) 
inside the SO coupled region. From the discrete version 
of the Schrodinger equation, one can show that the am- 
plitudes of the wave function in the SO coupled region 
and at those boundary lattice sites of the external leads 
connected directly to the SO coupled region will satisfy 
the following coupled equations: 

i?VT^(r3) = E H,ysirsa\vy)rT{v'^) 

- E h'^^r^-np'y'^T^i^^y'p')^ (5a) 

p 

-EV.^r..,v„,Vr'^(r.), (5b) 

where (rs, r(,) denote two nearest-neighbored lattice sites 
in the SO coupled region and (x^, ,Xp,) two nearest- 
neighbored boundary lattice sites in lead p' . ( For sim- 
plicity of notation, in the subscript of the Kronecker 
(5— function we have used simply a symbol Upiy, to de- 
note a boundary lattice site in the SO coupled region 
which is connected directly to a boundary lattice site 
Xp, — {l,y'p,) in lead p'. ) The matrix elements 

HsysijCsCT .v^a ) and Hp>{Xp,a ,Xp,<T ) can be written 
down directly from the Hamiltonian (1). Eq.(5a) and 
(5b) are the match conditions of the scattering wave func- 
tion on the borders between the SO coupled region and 
the leads, from which both the scattering wave function 
in the entire structure and all scattering amplitudes can 
be obtained simultaneously. Some details of the deriva- 
tions are given in the appendix. 

III. SOME RIGOROUS PROPERTIES OF 
EQUILIBRIUM STATES 

A controversial issue encountered in the study of spin- 
polarized transports in intrinsically SO coupled systems 
is that wether there can exist nonvanishing equilibrium 
background spin currents in such systems. Recently 
Rashba pointed out that, in a bulk two-dimensional elec- 
tron gas ( 2DEG ) with Rashba SO coupling, a finite equi- 
librium background spin current could be obtained if the 
conventional definition of spin current is applied to such 
systems ;S If this equilibrium background spin current 
does exist, it would imply that nonvanishing equilibrium 
spin polarizations should also exist near the edges of such 
a system due to the flow of the equilibrium background 
spin current. It was argued in Ref. [8] that such equilib- 
rium background spin currents are an artefact caused by 
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the improper use of the conventional definition of spin 
current to an intrinsically SO coupled system, i.e., the 
conventional definition of spin current cannot be applied 
in the presence of intrinsic SO coupling. Since it seemed 
that no consensus had been arrived on whether there is 
a uniquely correct definition for spin current in a SO 
coupled system2Si22ii2£i^, it would be meaningful if this 
controversial issue can be investigated in a somewhat dif- 
ferent way. In this section we will use the scattering 
wave function approach introduced in section II to inves- 
tigate wether there can exist nonvanishing equilibrium 
background spin currents and/or nonvanishing equilib- 
rium spin polarizations in a multi-terminal mesoscopic 
SO coupled system. Based on some simple but rigorous 
arguments, we will show that no finite equilibrium spin 
polarizations and/or finite equilibrium terminal spin cur- 
rents can exist in a multi-terminal mesoscopic SO coupled 
system. 

A. Absence of equilibrium spin polarizations 

In the tight-binding representation, the operator for 
the local spin density at a lattice site i reads 

k^ = \Y.^l^-Ap- (6) 

( For simplicity of notation, from now on we will use 
simply a symbol i to denote a lattice site in the en- 
tire structure, i.e., both in the SO coupled region and 
in the external leads. ) Under time reversal transfor- 
mation, the local spin density operator will transform as 

S{i) — > TS(i)T^^ = and the spin operator trans- 

form as aaf3 fSapT^^ = (-l)"+''(7*^ = whcrc 
a = —a, ^ = — /3, T = lOyK denotes the time-reversal 
transformation operator and K the conjugate operator. 

Within the framework of the standard Landauer- 
Biittiker's formalism, any physical quantities of a meso- 
scopic system are contributed to by all scattering states 
of conduction electrons incident from all contacts. These 
scattering states constitute an ensemble which can be 
specified by a chemical potential \Xp for each con- 
tact through their separate Fermi distribution function 
f{E,fip). For the problems discussed in the present pa- 
per, this ensemble consists of all scattering states de- 
scribed by the scattering wave functions given 
by Eqs.(3-4). In order that there is only one particle 
feeding into each incident channeP-^ when we use this 
ensemble of scattering wave functions to calculate the 
expectation value of an operator, one should first nor- 
malize the scattering wave function ipP"^'^ by a factor of 
l/VX ( L ^ oo is the length of lead p ), correspond- 
ing to that one changes the incident wave functions from 
gifc^a:p Q^Kr,^v I ^fl^_ By use of the ensemble of the 
normalized scattering wave functions {'0^""'^} and notic- 
ing that the density of states ( DOS ) for the mth trans- 
verse mode of lead p is given by ^ ^ = — , where 



Wpm = sin(fc'^) is the longitudinal velocity of the mth 
transverse mode of lead p, then one can see that the local 
spin density at a lattice site i ( either in the SO coupled 
region or in the external leads ) will be given by 

■pma ^ 

Q,/3 

where is the scattering wave function correspond- 

ing to an incident electron from the (ma) channel of lead 
p with a given energy E. This formula is valid both in 
the equilibrium and in the nonequilibrium states. If the 
system is in an equilibrium state, the chemical potential 
/ip will be independent of the lead label, i.e., = /i 
and f{E,iip) = f{E,^). Then in Eq.(7) the summation 
J2pma [• ■ ■] '^^'^ be performed first before carrying out the 
integration over energy E, and the result for this sum- 
mation can be expressed as 

E E [V'r^*(0(^'?)Q/3^r''(*)/H™ + H.c] 

afH pma 

Q/3 

= iY,{[G^{E)-G^{E)lp,^{^aU~H.C.}. (8) 

Q/3 

Here G^'^^iE) = [EI~H ± iO+]-^ is just the re- 
tarded and advanced Green's functions for the sys- 
tem, whose explicit spin-resolved matrix forms are given 

by GtU^) = i:p'm'a'ry'''*{^4"''"\j)/[E-E' ± 

iO^]^^ ( E is the incident energy correspond- 
ing to a scattering wave function ipP ™ ) ; and 

i[G^[E) — G^{E)\jp^ia is the spin-resolved spectral func- 
tion. If the total Hamiltonian H for the entire system 
is time- reversal invariant, the retarded and advanced 
Green's functions can be related by the time reversal 
transformation as 

Gt,p = {fG^f-^U,p = {-^r^'Gfl^-p (9) 

Combining Eq.® and Eq.® and taking into account 
the fact that fijapT'^ = (— 1)"+^(t*^ = — (Ta/3, one gets 
immediately that the right-hand side of Eq.® should 
vanish exactly, thus the spin density given by Eq. (7) van- 
ishes exactly in the equilibrium state, suggesting that no 
finite equilibrium spin polarizations can survive at any 
lattice site i in the entire structure ( both in the SO cou- 
pled region and in the leads ) . It should be noted that in 
arriving at this conclusion we have only made use of the 
assumption that the total Hamiltonian H for the entire 
system is time-reversal invariant ( which should be the 
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case in the absence of magnetic fields ) and did not in- 
volve the actual form of the SO coupling in the system, 
so it is a much general conclusion. 



B. Absence of equilibrium terminal spin currents 

In this subsection we discuss whether there can ex- 
ist nonvanishing equilibrium terminal spin currents in 
a multi-terminal mesoscopic SO coupled system. Since 
we have assumed that the leads are ideal and nonmag- 
netic ( i.e., described by a simple Hamiltonian Hp — 
—tpYj^-p.^-p.ycr{C'^.^Cpja + H.C.) ), the conventional defi- 
nitions of charge and spin currents can be well applied in 
the leads without ambiguities. According to the conven- 
tional definitions and in the lattice representation, the 
charge current and spin current ( with spin parallel to 
the a axia^ ) flowing from a lattice site Pi to a nearest- 
neighbored lattice site Pj in lead p can be given by the 
the corresponding particle density current as following. 



= -\J+ 



(10a) 



(10b) 



where Ip^p.^p. denotes the charge current operator and 
Ipp.^p. the spin current operator ( with spin parallel to 

the a axis ) and ^^p.^p the spin-resolved particle den- 
sity current operator and cr = ± denotes the spin-up 
and spin-down states with respect to the a axis. From 
the Heisenberg equation of motion for the on-site particle 
density: f^Np^ = ^[Np^,Hp], where iVp, = Cl^^Cp^a is 
the the on-site particle density operator in lead one can 
show easily that the spin-resolved particle density current 
flowing from a lattice site Pi to a nearest-neighbored lat- 
tice pj in lead p will be given by 



PjCr^PiO- 



^PiCr^PjO- 



)• (11) 



Now we calculate the terminal charge and spin cur- 
rents flowing along the longitudinal direction of a lead q. 
Firstly we consider the contribution of an incident elec- 
tron from the (mcr) channel of lead p to the longitudinal 
charge and spin currents ( with spin parallel to the a axis 
) flowing through a transverse cross-section ( saying, e.g., 
the cross-section at a; = a:^ ) of lead g, which by definition 
will be given by 



q/pmcr 



|^pmcr|2 



(12a) 



q Ipma 



{Xq -f l,yq) 



Vq 



h 

4^ 



(12b) 



where {xq,yq) and {xq + l,yq) denote two nearest- 
neighbored lattice sites along the longitudinal direction 
of lead q, -^lipP™"' {xq,yq)) denotes the normalized scat- 
tering wave function in lead q corresponding to the in- 
cident electron from the (mcr) channel of lead p ( which 
is given by Eqs.(3-4) ), and Vqn — -^sin(/c^) denotes 
the longitudinal velocity of the nth transverse mode in 
lead q and Vpm = -jf- sin(fc^) the longitudinal velocity 
of the mth transverse mode in lead p. The summation 
over the transverse coordinate yq runs over from 1 to Nq 
( Nq is the width of lead q )^, and the following or- 
thogonality relations for transverse modes in lead q have 
been applied in obtaining the last lines of Eq.(12a) and 
(12b): Y.y,x'L{yq)x'k{Vq) = Smn- K should be noted 
that, ii p = q, the results given by Eq.(12a) and (12b) 
will denote actually the contribution of an incident elec- 
tron from lead q to the charge and spin currents flowing 
in the same lead and { a' — ± ) denote actually 

the spin-flip { a a ) and non-spin-fiip { a = a ) 
reflection amplitudes. ( See the explanations given to 
Eqs.(3-4) in section II ). In such cases, the results given 
by Eq.(12a) and (12b) can be expressed as the subtrac- 
tion of the contributions due to the incident wave ( i.e., 
the terms proportional to Spq in Eq.(12a) and (12b) ) 
and the contributions due to the spin-flip and non-spin- 
flip reflected waves. 

The total terminal charge current Iq and the total ter- 
minal spin current Iq flowing in lead q will be obtained 
by summing the contributions of all incident electrons 
from all leads with the corresponding density of states ( 
see the explanations given above Eq.(7) ). Then we get 
that 

= tT. [dEf{E,f.p)Lj:\€Z'\'^-W 

J , t^om 

pma" n,a' ^ 

= JE/ dE]{E,p.p)\Tll,{E)-8pq8,,.Nq{E)-\ 

E / dE if{E, f,p)T^:, (E) - f{E, ^q)Tll. (i?)] , 



pacT 



(13a) 



= E /di?/(i?,Mp)Ei^(i<»+P-iC-n] 

= ^E/ dEfiE,^,p)[T^liE)^T^l{E)], (13b) 
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where r;;,(i?)^ ^ 



denotes ( by definition 



) the transmission probabiHty fi'om lead p with spin a to 
lead q with spin a' ( see Ref. [33] and also the explanations 
given in the appendix A ), and in obtaining the last line 
of Eq.(13a) the following relation has been applied^: 



(14) 



where Nq{E) is the total number of conducting trans- 
verse modes in lead q corresponding to a given energy 
E. As was discussed in detail in Ref. [33], this relation 
follows directly from the unitarity of the S-matrix, which 
is essential for the particle number conservation. 

Eq.(13a) is just the usual Landauer-Biittiker formula 
for terminal charge currents in a multi-terminal meso- 
scopic system'^'^. The second line in Eq.(13a) indicates 
clearly that the terminal charge current flowing in lead q 
can be expressed as the subtraction of the contributions 
due to all incident modes ( corresponding to the terms 
proportional to dpq ) and the contributions due to all out- 
going modes ( corresponding to the terms proportional 
to T^^, ) , which include both the transmitted waves from 
other leads { p ^ q ) and the reflected waves in lead q 
{ p = q ), noticing that T^^, denotes actually the spin- 
flip or non-spin-flip reflection probabilities from lead q 
to lead q ii p = q. Eq.(13b) is somewhat different from 
Eq.(13a) in appearance, but the terminal spin current 
given by Eq.(13b) can still be divided into two different 
kinds of contributions, namely the contributions due to 
the transmitted waves from other leads ( corresponding 
to those terms with p q hi the summation X^ptrl- ■ ■] ) 
and the contributions due to the spin-flip and non-spin- 
flip reflected waves in lead q ( corresponding to those 
terms with p — q in the summation X^po-I- ■ ■] )^ 

Eq.(13a) and (13b) are valid both in the equilibrium 
and in the nonequilibrium states. In the equilibrium 
state, since fip = ji { independent of the lead label ) 
and f{E,^p) = f{E,^), in Eq.(13a) and (13b) the sum- 
mation X^pcr [■ ■ ■] '■^^^ performed flrst before carrying 
out the energy integration and we get that 

Iq^^J dEf{E,^)J2^T^:'iE)-T^:,m, (ISa) 
I^^l- j dEf{E,^,)Y,[T^+{E)-Tll{E)]. (15b) 

pa 

Then by use of Eq.(14) one can see clearly that both ter- 
minal charge currents and terminal spin currents will van- 
ish exactly in the equilibrium state. It should be stressed 
that in arriving at this conclusion we did not involve 
the controversial issue of what is the correct deflnition of 
spin current in the central SO coupled region at all, so 
those ambiguities that might be caused by the use of an 
improper deflnition of spin current to the SO coupled re- 
gion have been eliminated in our derivations. Though we 



cannot prove that the spin current also vanishes exactly 
inside the SO coupled region based on the approach ap- 
plied above, however, for a mesoscopic system only the 
terminal ( charge or spin ) currents are the real useful 
quantities from the practical point of view ( i.e., one 
need to add external contacts to induct the charge or 
spin currents out of a mesoscopic sample ). We note 
that a similar conclusion as was obtained above has also 
been derived in Ref. [37] based on some somewhat differ- 
ent arguments. Compared with the derivations given in 
Ref. [37] , the arguments given above seem to be more sim- 
ple and more transparent in principle. It also should be 
noted that, based on a similar Landauer-Biittiker formal- 
ism, it was argued in Ref. [38] that the equilibrium termi- 
nal spin currents should indeed take place in a three ter- 
minal system with spin-orbit coupling^^, in contradiction 
to the conclusion obtained in the present paper and in 
Ref. [37]. To our understandings, this contradiction was 
caused by the fact that the contributions due to the spin- 
flip and non-spin-flip reflections in the leads ( induced by 
the scatterings from the central SO coupled region ) was 
neglected in the calculations of terminal spin currents 
performed in Ref. [38]. In contrast, in the calculations of 
terminal spin currents performed in the present paper, 
the contributions due to the spin-flip and non-spin-flip 
reflections in the leads induced by the scatterings from 
the central SO coupled region have been treated in an ac- 
curate and strict way, assuming that the leads are ideal 
and nonmagnetic. 



IV. SOME SYMMETRY PROPERTIES OF 
NONEQUILIBRIUM SPIN CURRENTS AND 
SPIN POLARIZATIONS IN TWO-TERMINAL 
MESOSCOPIC SO COUPLED SYSTEMS 



When a multi-terminal mesoscopic SO coupled system 
is driven into a nonequilibrium state ( i.e., there is charge 
current flow between different leads ), nonequilibrium 
spin polarizations and/or terminal spin currents may be 
induced by the charge current flow. In this section we 
discuss some symmetry properties of such nonequilibrium 
spin currents and spin polarizations. For clarity, we take 
a typical two-terminal mesoscopic structure as shown in 
Fig.?? as the example, where a ballistic two-dimensional 
electron gas ( 2DEG ) with Rashba and/or k- linear Dres- 
selhaus SO coupling is attached to two ideal leads. 



lead 1 










lead 2 




O X 




2DEG 





FIG. 2: Schematic geometry of a two-terminal mesoscopic SO 
coupled system. 
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A. Non-antisymmetric lateral edge spin 
accumulations 

The study of nonequilibrium lateral edge spin accu- 
mulation induced by a longitudinal charge current in a 
thin strip of a two-dimensional electron gas with intrinsic 
SO coupling is of great theoretical interest because of its 
close relations with the intrinsic spin Hall effect in such 
systems. It was generally believed that the principal ob- 
servable signature of the intrinsic spin Hall effect in a SO 
coupled system is that, when a longitudinal charge cur- 
rent circulates through such a system with a thin strip 
geometry, antisymmetric lateral edge spin accumulation 
( polarized perpendicular to the 2DEG plane ) will be 
induced at the two lateral edges of the strip due to the 
flow of the transverse spin Hall current. Several recent 
numerical calculations had demonstrated that a longitu- 
dinal charge current circulating through a thin strip of 
a ballistic two-dimensional electron gas with Rashba SO 
coupling does can lead to the generation of antisymmet- 
ric edge spin accumulation at the two lateral edges of the 
strip, and the antisymmetric character of the transverse 
spatial distribution of the lateral edge spin accumulation 
( i.e., {Sz{x,y)) = —{Sz{x, —y)) ) had been argued to be 
a strong support of the existence of intrinsic spin Hall 
effect in such mesoscopic SO coupled systems^-. Here 
we discuss this issue from a different point of view. We 
will show that, when a longitudinal charge current circu- 
lates through a thin strip of a ballistic two-dimensional 
electron gas with both Rashba and Dresselhaus SO cou- 
pling, the transverse spatial distribution of the induced 
nonequilibrium lateral edge spin accumulation ( polar- 
ized perpendicular to the 2DEG plane ) are in general 
non-antisymmetric. The non-antisymmetric character 
of the lateral edge spin accumulation contradicts seri- 
ously with the usual physical pictures of spin Hall effect, 
though according to some theoretical predictions, intrin- 
sic spin Hall effect should also survive in the presence of 
both Rashba and Dresselhaus SO couplingi^. The non- 
antisymmetric character of the lateral edge spin accumu- 
lation implies that, in addition to the intrinsic spin Hall 
effect, there may exist some other physical reasons that 
might also lead to the generation of nonequilibrium lat- 
eral edge spin accumulation in a SO coupled system ( 
with a thin strip geometry ) when a longitudinal charge 
current circulates through iti^i^i 

Firstly let us look at what symmetry relations can be 
obtained for the nonequilibrium lateral edge spin accu- 
mulation induced by a longitudinal charge current based 
on the symmetry analysis of the Hamiltonian of the sys- 
tem under study ( sketched in Fig. 2 ). If only Rashba ( 
or only Dresselhaus ) SO coupling presents, based on the 
symmetry properties of the Hamiltonian of the structure 
under study, one can show rigorous that the nonequi- 
librium lateral edge spin accumulation does should be 
antisymmetric at the two lateral edges. Let us consider 
first the case in which only Rashba SO coupling presents 
( i.e., the Dresselhaus SO coupling strength is zero ). 



If only Rashba SO coupling presents, from Eq.(2a) and 
Fig. 2 one can see that the Hamiltonian of the entire 
structure is invariant under the combined transforma- 
tion of the real space reflection y ^ —y and the spin 
space rotation around the Sy axis ( with an angle tt 
). From this invariance one can get immediately that 
(Sz{x,y))i = -{Sz{x,-y))i, where {Sz)i denotes the 
nonequilibrium spin density induced by a longitudinal 
charge current flowing from lead 1 to lead 2. It is inter- 
esting to note that this antisymmetric relation can be de- 
duced directly from the symmetry of the structure under 
study but without need to resort to the concept of spin 
Hall effect at all. Next, let us consider the case in which 
only Dresselhaus SO coupling presents ( i.e., the Rashba 
SO coupling strength is zero ). If only Dresselhaus SO 
coupling presents, then from Eq.(2b) and Fig. 2 one can 
see that the Hamiltonian of the entire structure is invari- 
ant under the combined transformation of the real space 
reflection y ^ ~y and the spin space rotation around the 
Sx axis ( with an angle tt ). From this invariance one also 
gets immediately that {Sz{x,y))i = —{Sz{x,—y))i, i.e., 
the nonequilibrium lateral edge spin accumulation still 
should be antisymmetric at the two lateral edges if only 
Dresselhaus SO coupling presents. 

If both Rashba and Dresselhaus SO couplings are 
present, then from Eqs.(2a-2b) and Fig. 2 one can see that 
the total Hamiltonian of the entire structure is invariant 
under the combined transformation of the real space cen- 
ter inversion r =^ — r and the spin space rotation around 
the Sz axis ( with an angle tt ). From this invariance 
one can get that {Sz(x,y))i = {Sz{—x,—y))^i, where 
{Sz)-i denotes the nonequilibrium spin density induced 
by a longitudinal charge current flowing from lead 2 to 
lead 1. On the other hand, from Eq.([7]) one can see that 
in the linear response regime one has 

{Sz{x,y))i = ~{Sz{x,y))-i. (16) 

Combining the two relations {Sz{x,y))i = 
{Sz(-x,-y))^i and {Sz{x,y))i = -{Sz{x,y))-i , 
then the following symmetry relation can be obtained 
for the nonequilibrium spin accumulation induced by 
a longitudinal charge current flowing from lead 1 to 
lead 2: {Sz{x,y))i = -~{Sz{-x,-y))i. This symmetry 
relation implies that the transverse spatial distribution 
of the nonequilibrium lateral edge spin accumulation 
will be antisymmetric in the center cross-section of the 
strip, i.e., (^^(O,?/))/ = -(5z(0, -y))/. Due to the 
existence of this symmetry relation, one can deduce that 
in an infinite strip ( i.e, the length of the strip tends to 
infinity and hence the effects of the contacted leads can 
be neglected ), the transverse spatial distribution of the 
nonequilibrium lateral edge spin accumulation will still 
be antisymmetric ( i.e., {Sz{x,y))i — —{Sz{x,—y))i for 
all X ) in the presence of both Rashba and Dresselhaus 
SO coupling. However, unlike the case in which only 
Rashba ( or only Dresselhaus ) SO couphng presents, 
for a thin strip of finite length, in the presence of both 
Rashba and Dresselhaus SO coupling, one cannot de- 
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duce a general antisymmetric relation for the transverse 
spatial distribution of the nonequilibrium lateral edge 
spin accumulation based on the symmetry analysis of 
the total Hamiltonian of the entire structure under 
study. From the theoretical points of view, this is 
due to the fact that, in the presence of both Rashba 
and Dresselhaus SO coupling, the total Hamiltonian of 
the entire structure under study ( i.e., a thin strip of 
finite length contacted to two ideal leads ) is no longer 
invariant under the combined transformation of the real 
space reflection y =^ —y and the spin space rotation 
around the Sy { or Sx ) axis with an angle tt. Indeed, as 
will be shown below, this non-antisymmetric character 
can be verified by detailed numerical calculations. 

One particular interesting case is that the Rashba and 
the Dresselhaus SO couphng strengths are equal ( i.e., 
tR — tu or tn — —to ). In this particular case, the 
total Hamiltonian is invariant under the following uni- 
tary transformation in spin space ( while the real space 
coordinate r remain unchanged ): U^HU+ = H { ii 
tR = to ) or JJ^HJJ- = H { ii tR = -to ), where 
?7+ = {ax + ay) 1^/2 and [/_ — [ax —ay)/^/2. Under this 
unitary transformation, the spin operators will transform 
as following: a^ — > -ct^, cFx ^ ay { if Ir = to ) or 
ax ^ —ay ( \i tR = —tn ). Since the real space coordi- 
nate r remain unchanged under this symmetry manipu- 
lation, from the above symmetry properties in spin space 
one gets immediately that {Sz{x,y))i — —{Sz{x,y))i, 
suggesting that {Sz{x,y))i should vanish everywhere if 
the Rashba and the Dresselhaus SO coupling strengths 
are equal. This conclusion is in exact agreement with 
the corresponding numerical results obtained based on 
the scattering wave approach introduced in section II- 
ni, which shows that {Sz{x,y))i does vanish everywhere 
in the particular case of tR = to or tR = —to- 

To show more explicitly the non-antisymmetric charac- 
ter of the lateral edge spin accumulation induced by a lon- 
gitudinal charge current in a 2DEG strip of finite length 
with both Rashba and Dresselhaus SO coupling, in Fig. 3 
we plotted a typical pattern of the two-dimensional spa- 
tial distribution of the nonequilibrium spin density (Sz) 
in the strip obtained by numerical calculations with the 
scattering wave function approach introduced in Sec. II- 
III. In our numerical calculations we take the typical val- 
ues of the electron effective mass m = 0.04me, the lattice 
constant a = 3nm, and the 2DEG strip contains 120 x 40 
lattice sizes. The chemical potentials in the two leads 
are set by fixing the longitudinal charge current to flow 
from lead 1 to lead 2 as shown in Fig. 2 and fixing the 
longitudinal charge current density to 100 fi A /1. 5 fj,m ( as 
reported in Ref.[26] ). From Fig. 3 one can see clearly that 
the transverse spatial distribution of the nonequilibrium 
spin density (Sz) in the strip is non-antisymmetric in gen- 
eral ( i.e., {Szix, y))i ^ -{Sz{x, -y))i for general a; ), ex- 
cept in the center cross-section ( i.e., a; = ) of the strip. 
The non-antisymmetric character of the transverse spa- 
tial distribution of the nonequilibrium spin density {Sz) 
can be more clearly seen from Fig. 4(a), where we plotted 



(S^lx.y)), (h/4x) 
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FIG. 3: (Color online) A typical pattern of the two- 
dimensional spatial distribution of the current induced 
nonequilibrium spin density {Sz) in a two-terminal structure ( 
sketched in Fig. 2 ) in the presence of both Rashba and Dressel- 
haus SO coupling. The Rashba and Dresselhaus SO coupling 
strength is set to tn/t = 0.08 and tn/t = 0.02. 



several typical patterns of the profiles of the transverse 
spatial distributions of the nonequilibrium spin density 
{Sz) in a cross-section of the strip at x ^ 0. ( For 
comparison, the corresponding results obtained in the 
case that only Rashba or only Dresselhaus SO coupling 
presents were also plotted in Fig. 4(b) ). The three typi- 
cal patterns shown in Fig. 4(a) are obtained by fixing the 
Dresselhaus SO coupling strength to to — 0.02t ( t is the 
spin-independent hopping parameter ) and varying the 
Rashba SO coupling strength tR. From Fig. 4(a) one can 
see clearly that, the transverse spatial distribution of the 
nonequilibrium spin density {Sz{x,y))i can have either 
the same signs or opposite signs at the two lateral edges 
of the strip, depending on the ratios of tR/to- Even in the 
case that {Sz) has opposite signs at the two lateral edges, 
the transverse spatial distributions of {Sz) may still not 
be antisymmetric ( i.e., {Sz{x, y)) ^ —{Sz{x, —y)) ), con- 
tradicting significantly with the usual physical pictures of 
spin Hall effect. The non-antisymmetric character of the 
lateral edge spin accumulation suggests that some cau- 
tions may need to be taken when attributing the nonequi- 
librium lateral edge spin accumulation induced by a lon- 
gitudinal charge current in a thin strip of a SO coupled 
system to a spin Hall effect, especially in the mesoscopic 
regime. 

The results shown in Figs. 3-4 are obtained in the ab- 
sence of impurity scatterings. One can show that the 
symmetry properties shown in Fig. 3-4 are robust against 
spinless weak impurity scatterings. To model spinless 
weak disorder scatterings, we assume that the on-site 
energy at lattice sites in the 2DEG strip are randomly 
distributed in a narrow energy region [— VF, W^], where 
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FIG. 4: (Color online) (a) Some typical profiles of the trans- 
verse spatial distributions of the nonequilibrium spin density 
(Sz) induced by a longitudinal charge current in a 2DEG strip 
with both Rashba and Dresselhaus SO couplings, which shows 
clearly that the transverse spatial distributions of (Sz) are 
non-antisymmetric in general at both edges of the strip in 
the presence of both Rashba and Dresselhaus SO coupling. 
{Sz} vanishes everywhere in the particular case oi ta = to or 
tn = —to { not shown explicitly in the figure ). (b) The corre- 
sponding results obtained in the case that only Rashba or only 
Dresselhaus SO coupling presents, which shows clearly that 
the transverse spatial distributions of (Sz) are antisymmetric 
at both edges of the strip if only Rashba or only Dresselhaus 
SO coupling presents. 



W is the amplitude of the on-site energy fluctuations 
characterizing the disorder strength. ( In the absence 
of disorder scatterings, the on-site energy at each lattice 
site can be set simply to zero. ) We calculate the spin 
density for a number of random impurity configurations 
and then do impurity average. In Fig. 5 we show the 
variations of the transverse spatial distributions of the 
nonequilibrium spin density (Sz) in a cross-section of the 
strip as the disorder strength increases, from which one 
can see that the symmetry properties of the transverse 
spatial distributions of the nonequilibrium spin density 
are robust against spinless weak disorder scatterings. 



o.ooe 

0.006 

5 0.004- 

a; 0.002 

° 0.000 

5^ -0-002 
V 

-0.004 
-0.006 
-0.008 



t /t=0.08, t„/t=0.02 




FIG. 5: (Color online)The profiles of the transverse spatial 
distributions of the nonequilibrium spin density {Sz} in the 
presence of disorder. We have done impurity average over 
1000 random impurity configurations for each case. 



B. Non-conservative terminal spin currents 

As has been mentioned in the introduction section, a 
much controversial issue related to the study of spin- 
polarized transports in SO coupled systems is that 
whether spin currents should be a conserved quantity 
and what is the correct definition of spin currents in such 
systems. It was now well established that, in the pres- 
ence of SO coupling, spin current calculated based on 
the conventional definition is not a conserved quantity, 
and in order to make spin current a conserved quantity 
in the presence of SO coupling, significant modifications 
to the conventional definition will be neededi^i^S Never- 
theless, it seemed that no consensus had been arrived on 
whether spin current should be a conserved quantity in 
a SO coupled system or whether there is a uniquely cor- 
rect definition for spin current in such a system.— i^i^!^ 
Below we will discuss this controversial issue from a dif- 
ferent point of view, i.e., we do not consider the problem 
that what is the correct definition of spin current in a 
SO coupled system but focus our discussion on the ques- 
tion that whether the terminal spin currents in a multi- 
terminal mesoscopic SO coupled system are conservative. 
As mentioned earlier, for a mesoscopic SO coupled sys- 
tem, only the terminal spin currents are the real useful 
quantities. By use of the two-probe mesoscopic struc- 
ture shown in Fig. 2 as the example, we will show explic- 
itly that the terminal spin currents in a multi-terminal 
mesoscopic SO coupled system are non-conservative in 
general, i.e., the total spin currents flowing into the SO 
coupled region are not equal to the total spin currents 
flowing out of the same region. To illustrate this point 
clearly, we take a two-terminal mesoscopic system with 
both Rashba and Dresselhaus SO coupling as the exam- 
ple. In Fig. 6 (a) and (b) we plotted the terminal spin 
currents If and /| in the two leads ( with spin parallel 
to the z axis ) and the terminal spin currents If and in 
the two leads ( with spin parallel to the y axis ) as a func- 
tion of the Rashba SO coupling strength i^, respectively. 
In our calculations we fix the Dresselhaus SO coupling 
strength to = 0.02t and fix the longitudinal charge 
current density to 100/xA/1.5/im. The lattice constant 
a = 3nm and the lattice size of the 2DEG strip is taken 
to be 100 X 40. The positive direction of the spin current 
flow is defined to be from lead 1 to lead 2. From Fig. 6(a) 
one can see that the terminal spin currents If and /| in 
the two leads have the same signs, which means that the 
terminal spin currents with spin parallel to the z axis will 
flow from lead 1 into the SO coupled region and then flow 
out of the SO coupled region into lead 2, similar to the 
usual charge current transport. From Fig. 6(b) one can 
see that the terminal spin currents If and /| in the two 
leads have opposite signs, which means that the terminal 
spin currents with spin parallel to the y axis will flow 
out of the SO coupled region in both leads and hence 
are non-conservative ( i.e., the spin current flowing into 
the SO coupled region does not equal to the spin current 
flowing out of the same region ) . Similarly one can show 
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FIG. 6: (Color online) (a) The terminal spin currents and 
— /f ( divided by the voltage ) as a function of the Rashba SO 
coupling strenght tu ( in units of t ). (b) The terminal spin 
currents If and /| ( divided by the voltage ) as a function of 
the Rashba SO coupling strenght tn. The figures show that 
the terminal spin currents Ii and /| in the two leads have the 
same signs and the terminal spin currents If and If in the 
two leads have opposite signs. ( Note that for clarity a minus 
sign is added before If in Fig. 6(a) ). The parameters used 
are given in the text or shown in the figures. 



that the terminal spin currents with spin parallel to the 
X axis have also opposite signs in the two leads, similar 
to the case shown in Fig. 6(b). This simple example il- 
lustrates explicitly that the terminal spin currents in a 
multi-terminal mesoscopic SO coupled system are non- 
conservative in general. It should be stressed that this 
non-conservation of terminal spin currents is not caused 
by the use of an improper definition of spin current but is 
intrinsic to spin-dependent transports in mesoscopic SO 
coupled systems. In fact, in our calculations we did not 
involve the controversial issue of what is the correct def- 
inition of spin current in the SO coupled region at all, so 
the ambiguities that may be caused by the use of an im- 
proper definition of spin current to the SO coupled region 
have been eliminated in our calculations. 



CONCLUSION 



In summary, based on a scattering wave function ap- 
proach, in this paper we have studied theoretically some 
symmetry properties of spin currents and spin polariza- 
tions in a multi-terminal mesoscopic structure in which 
a spin-orbit coupled system is contacted to several ideal 
and nonmagnetic external leads. Some interesting new 
results were obtained based on the symmetry analysis 
of spin currents and spin polarizations in such a multi- 
terminal mesoscopic structure. First, we showed that in 



the equilibrium state no finite spin polarizations can ex- 
ist both in the leads and in the central SO coupled region 
and also no finite equilibrium terminal spin currents can 
survive. Second, we showed that the lateral edge spin ac- 
cumulation induced by a longitudinal charge current in a 
thin strip of a ballistic two-dimensional electron gas with 
both Rashba and Dresselhaus SO coupling may be non- 
antisymmetric in general, which implies that some cau- 
tions may need to be taken when attributing the nonequi- 
librium lateral edge spin accumulation induced by a lon- 
gitudinal charge current in a thin strip of such a system 
to a spin Hall effect, especially in the mesoscopic regime. 
Finally, by use of a typical two-probe structure as the 
example, we showed explicitly that the nonequilibrium 
terminal spin currents in a multi-terminal mesoscopic 
SO coupled system may be non-conservative in general. 
Some symmetry properties discussed in the present pa- 
per might also be helpful for clarifying some controversial 
issues related to the study of spin-dependent transports 
in macroscopic SO coupled systems. 
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APPENDIX A: SOME DETAILS FOR THE 

DERIVATIONS OF THE SCATTERING 
AMPLITUDES AND THE TRANSMISSION 
PROBABILITIES 

In this appendix we give some details on how to derive 
the scattering amplitudes and the transmission probabil- 
ities from the scattering wave function approach intro- 
duced in Sec. II. For convenience of notation, we arrange 
the scattering wave function ipP™'^{ri) inside the SO cou- 
pled region into a column vector whose dimension is 
2N { N is the total number of lattice sites in the SO 
coupled region ) and arrange the scattering amplitudes 
'f'qna' ^ column vcctor $ whose dimension is 2M ( 
M — J2p^p ^iiti is the width of lead p ). Substi- 
tuting Eqs.(3-4) into Eqs.(5a-5b) and making use of the 
orthogonality relations for the transverse modes in the 
leads, one can show that the two column vectors and 
$ will satisfy the following relations: 



A*, 



b + B$,C$ = d + D^-, 



(Al) 



where A, B, C, D are four rectangular matrices with the 
dimensions of 2N x 2iV, 2N x 2M, 2M x 2Af, and 
2M X 2N, respectively; b and d are two column vec- 
tors with the dimensions of 2N and 2M, respectively. 
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The elements of these matrices and cohimn vectors can 
be written down expUcitly as 



A 

'B{np"ycr" ,p'ni'(j') 

'D{p'm'a' , npiiyo") 
Cip'm'a',p"m"a") 

d{p'm'a') 



= EI-H. 



sys, 



— Spii ^pidu" .a'tp' sXm' (j/p" ) J 
p' 5(j" .a' ^m" m'ip' •> 

<5pp' Saa'tps Xm (j/p' ) ^ 



6pp' ^ra'm' ^aa' 



(A2) 



where I stands for the identity matrix. The indices for 
leads, transverse modes, lattice sites and spins can take 
all possible values. ( For simplicity of notation, we have 
used simply a symbol Upiyi to denote a boundary lattice 
site in the SO coupled region which is connected directly 
to a boundary lattice site x^, = (1,?/^/) in lead p' . ) 
Eq. (|Aip is just a compact form of the match conditions 
(5a-5b) on the borders between the leads and the SO cou- 
pled region, from which both the scattering amplitudes 
^^"^ transmission probabihties {T^",} can be 
obtained readily. 

To derive a compact formula for the transmission prob- 
abilities between two leads, we define an auxiliary matrix 
= BC"^D. By use of Eq. (U2l) the matrix elements 
of E^ can be written down readily as following, 

^^inp,y,a\np.y^y) ^-Y^fl^P;^,(y,)^i,(y,)e<', 

ra ^ 

(A3) 

and all other matrix elements not shown explicitly above 
are zero. With the help of this auxiliary matrix, from 
Eg. ljAip one can get that 

= (A- E'^)"i(b + BC"M) = G-^g, (A4) 

where g is a column vector defined by g = b + BC^^d 
and is a matrix defined by = (A — E^)^^ = 



[El — Hgys — E-'^]"^, which is just the usual retarded 
Green's function. By use of Eq. (jA2p . the elements of 
the column vector g can also be written down readily as 
following. 



g(np'ycr') = 2iSpp,Saa'tpsSm(kP^)x^{y)- 



(A5) 



By substituting Eq. (|A4p into Eq. (|A1|) one gets that $ = 
C^^d + C^^DG-'^g. Inserting Eq.(A5) into this formula 
and making use of Eq.(A2), one can show readily that 
the scattering amplitudes <t)^l^Ji will be given by 



ipma- 
^p'ln' a' 



Spp' Sjrim' ^aa' ~^ ^^^p' tp^ s ^ ^ ^ps sin(/c^) 
^Vp' ' ^Vp )Xm{yp)- (A6) 



The total transmission probability of a conduction elec- 
tron from lead p ( with spin index a ) to lead p' ( with 



ipma 
^p'jn'a' 



spin index ct ) is defined by T^/'^/ = 

^ m.m' 

where Wp'm' — ji^^p' sin(fc^/) is the longitudinal velocity 
of the mode m' in lead p' . Substituting Eq. (|A6|) into this 
formula, for p ^ p' one can get that 



(A7) 



where G^^, and G^^,(= G^}^) are the spin-resolved re- 
tarded and advanced Green's functions, respectively, and 
^^iVp^Vp) is defined by 



^^{vp^Vp) = ^ij^)'^Xm{yp)vpm,xm{yp)- 



(A8) 



The transmission probabilities given by Eq. (|A7[1 have ex- 
actly the same form as was obtained by the usual Green's 
function approach^. 
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